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ABSTRACT 

We discover that some unstable vacua have long memory. By that we mean that even in the 
theories containing only massive particles, there are correllators and expectation values which grow 
with time. We examine the cases of instabilities caused by the constant electric fields, expanding 
and contracting universes and, most importantly, the global de Sitter space. In the last case the 
interaction leads to a remarkable UV/IR mixing and to a large back reaction. This gives reasons 
to believe that the cosmological constant problem could be resolved by the infrared physics. 



1 Introduction 

The problem of cosmological constant presents a serious challenge to modern physics. Recently 
the following physical mechanism was proposed to overcome this problem \2\ . Let us imagine that 
the bare cosmological constant is present in the lagrangian. As is well known, it is the cause of 
gravitational repulsion, resulting in the accelerated expansion of the universe (the contraction is 
also possible, but we will discuss it later). According to the proposal, this acceleration leads to the 
explosive particle production. The gravitational attraction between the particles slows down the 
acceleration and thus reduces (asymptotically to zero) the effective cosmological constant. 

There are many puzzles associated with this proposal. Particles in the curved space are ill 
defined, does it make sense to ascribe to them a real physical effect? Even if it does, the Universe 
is exponentially expanding, so it may seem that the particles get diluted; isn't their back reaction 
negligible? How can massive particles considered in [2^ lead to large infrared effects? 

In this paper we will try to provide some clarifications, as well as present some new results. It 
is helpful to consider various cases of unstable vacua and to make their comparative studies. The 
basic origin of the difficulties lies in the non-equilibrium quantum field theory and are common to 
all the cases. Let us demonstrate this with the following example. 

Consider a nucleus with a charge Z and hit it with a 7-quantum of energy lo which produces a 
pair e^e~ . The electron forms a bound state while the positron escapes to infinity. The threshold 
singularity in w is located at 

uj = 2m- \Eb\ 

where |£'b(Z)| is the binding energy. We see that when the nucleus becomes supercharged, 
\EB{Z)\ = 2m, we get a long-ranged correlations in time since the threshold is now at a;= 0. 



Below we will find that such "long- memory" is a crucial factor in the non-equilibrium dynamics. 
In field theory it leads to the "adiabatic catastrophe" [IJ and to the obstruction to Wick's rotation. 
In the next section we shall briefly summarize the situation. 

Another puzzle mentioned above is the dilution of particles in the expanding universe, the size 
of which grows as a{t) ^ e*. However, in any physical quantity, this effect always cancels with 
the exponentially growing number of the comoving modes. The covariant cut-off for the comoving 
momentum k is given by fc < /cmax MpianckCi{t) , and this causes the growth. The naive reason 
for the above compensation is that the change a — > Aa is a coordinate transformation and a 
dependence can't be physical. As will be discussed below, there are caveats to this argument, but 
on a qualitative level they are unimportant. 

What about the infrared effects generated by the massive particles? They are not related to 
the interaction of these particles, which is short-ranged as usual. Their origin lies in the fact 
that the original vacuum is unstable with the non-zero decay rate. Therefore we get perturbative 
corrections containing secular terms, which represent the fact that, as the time goes, it is less and 
less likely for the vacuum to remain intact. 



2 Lorentzian vs Euclidean calculations 

A crude example of these phenomena is provided by a hot plasma in a box. One can calcu- 
late various correlations using the Euclidean approach by introducing Matsubara's periodic time. 
The resulting physical quantities will describe a well-defined physics of plasma kept at fixed tem- 
perature by the external sources. However, if such sources are absent, the plasma will cool off. 
In this case the "Euclidean" calculation is inadequate and we have to use the Schwinger-Keldysh 
approach. 

Similar situation exists in the dS space. This space is defined by the equation 

^2 2 1 

n ~ Hq — 1 

It is tempting to define the following rules of the game. Make a "Wick rotation", no => ino, which 
transforms the dS space into a sphere. On a sphere massive particles will never generate any 
IR divergences, at least perturbatively. Suppose now that we calculate the correlation functions 
<(p(ni)... 93(71 Ar)> on a sphere and then continue them back to the dS space. Such a strategy was 
briefiy considered and rejected in [2J, but taken as a fundamental definition of the theory in [12j. 
Moreover, it was shown in these two papers that the analytic continuation from a sphere doesn't 
lead to any pathologies. 

Our claim is that, just as in the above case of plasma, this "Euclidean" approach describes 
the de Sitter space artificially kept at fixed Gibbons-Hawking temperature. This can be achieved 
by surrounding the dS patch by reflecting walls sending all emitted radiation back. However 
without these magic devices Euclidean approach is inadequate. 

A very similar situation occurs in the case of the Schwarzschild black holes — their Euclidean 
geometry is the well-known cigar and the field theory on the cigar describes eternal black hole in 
which the Hawking radiation is being returned to keep the equilibrium. Let us notice also that in 
the case of the black hole the Euclidean approach uses r > 2M region where (after Wick's rotation) 
the killing vectors are positive. 

As is well known, its Euclidean metric 

^ _ 2M ^ — ^4fT + ^'dn 
r ' (1-^) 

is complete and non-singular if f ^ t+8nM. It has a geometry of a cigar. According to Gibbons and 
Hawking, periodicity in imaginary time is the indication of Hawking's temperature. An interacting 
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field theory on this manifold is also well defined. We can get a set of Green's functions by an analytic 
continuation to real time. 

However, as well known, these continuation will describe the black hole on "life support" - 
its temperature must be kept constant. The real black hole evaporates and can't be described 
by Wick's rotation. The same is true for the dS space. A typical puzzle here is that while the 
black hole evaporates into the outer space, the dS space has nowhere to go. In fact, this puzzle is 
psychological. The dS space simply creates an avalanche of particles within itself . One can also 
visualize this by immersing a large patch of the dS space in the Minkowski space-time. In this case 
the created particles will populate the ambient space. For the proper description of the dS space 
one must use the Schwinger-Keldysh approach with real time. 

In the case of the Bunch-Davies vacuum there are some further puzzles. The wave functional 
in this case is of the Hartle-Hawking type: 

ip\aM=fixed 

where we integrate over fields on a hemisphere. The correlators are given by 

G=(V'o, (p{ni),...,ip{nN)ilJo)- 

The dS vector n can be parameterized as n = (sinht. (coshf)ra), the southern hemisphere is 
described by — i7r/2 ^ t ^ 0, while the northern one corresponds to i7r/2 ^ f ^ 0. If we consider a 
perturbation theory analytically continued to dS we have formally 

G{Euclidean) = (V'O, 7c'y?(ni)...(^(njv) 6"^ ^Vo) 

where the contour C for the correllators on a sphere is shown at Fig.l left, while the contour 
after analytic continuation to dS space is shown on the right. These contours must be used 



XK 
2 



2 




Figure 1 : The contours of integration on a sphere and after analytic continuation to dS space. 

in the Schwinger-Keldysh diagrams. One might think that this should give the same result as 
the standard Schwinger-Keldysh contour for the Lorentzian approach, Fig. 2. 

This is not the case. The reason for the discrepancy lies in the adiabatic catastrophe. We have 
to introduce the adiabatic switching of the interaction and let the system develop freely: 

G^Lorent.ian)ini,...) = (0|T^(ni)...(^(n;v)e-' ^(^""^'^''^"10) 

where £ — >■ 0. 

As we will show below the IR cut-off e does not disappear and that implies the instability 
of the dS space. The reason for the breakdown of adiabaticity is that there are states with negative 
energies and so the arbitrary slow external field can trigger pair production. 
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Figure 2: Standard contour of integration in Lorentzian approach. 

3 Electric fields and new anomaly 

Pair production by electric fields has been discussed in the hundreds of valuable papers, see ref- 
erences in [3j. We return to this problem for two reasons. First, we need to present it in the 
form, which can be easily generalized to the gravitational case. Second, we will find an unusual 
anomalous vacuum polarization which may have unexpected applications. 

Let us consider massive scalar field in electric field, described by a time-dependent vector 
potential Ai{t). We assume that electric field is switched on and off adiabatically. This means 
that it has the form E = E{-^) so that for \t\ << T, it remains constant while for \t\ >> T, 
E ^ 0. A good concrete example of such behavior (already considered in [4 ) is to take 



Ai{t) = ETtanh 



^(^) = — rj-rT2 J but explicit shape of the potential is not important. The Klein-Gordon equation 
has the form 

a) 



(a^ + {k^A{t)f + kl + m^)^ = 



We are interested in the 'in' solution which is defined as the 'Jost function', i.e. it has the asymptotic 
behavior 

1 



(Pin(t,k) 



2uj: 



(2) 



where — \J (k — A(±oo)) + fc^ + m?. The solution is normalized by the condition that 
Wronskian = 1. 

As we go to the late time t oo, we have 



Vin{t., k) -^f_ 



a{k)e 



me'' 



(3) 



where a and l3 are Bogolyubov coefficients also related to the transmission and refiection ampli- 
tudes. 

If we start with ipi„ and blindly apply the WKB approximation, we get 



^in{t, k) 



1 



-i / uik{t')dt' 

=e 



(4) 



for late times, with uJk{t) = \J [k — A{t)y' -I- fc^ -I- m^. Of course in this way we loose the over 
barrier reflection and thus the above formula can't be valid everywhere. Indeed the WKB requires 



4 



that the de BrogUe wave length A = — satisfies 



_ dX _ 
^ ^ ~dt ^ 



{k-A)A{t) 



< 1 



We see that the approximation is good for the early times when j/c — ^ m. In this case 



E 



E 



\k-A\'^ 



< 1 



However, around the point where the mode 'reaches the horizon', defined by fc = A{tk)^ we get 
7^1 and WKB breaks down. As we go to t ^ i^, j/c — starts growing again and the WKB is 
valid again. In this region it must contain two exponentials: 
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a{k)e ° + I3{k)e ° 



(5) 



As usual, a and /? can be found by matching ^ and ([5]). 

In the domain \t\ ^ T the electric field is constant and A{t) ~ Et. The equation ([T]) now 
depends on the variable t — hence ipm ~ fin{t ~ "§)• The function /„, as well known, is the 
parabolic cylinder function 



D_ 



E' 



t 



but we will not need their explicit form. What is important is that due to the symmetry fc — /c + k, 
i — > t — the resulting a and j3 do not depend on fc in a certain range, which we determine in a 
moment (but do depend on k± and m). 

To find this range we notice that the 'horizon crossing' (fc = A{t)) occurs at = -g. We can 
use the constant field approximation only if tk ^ T. Hence we conclude that a and /3 do not 
depend on k only if A{—oo) < k < A{oo). Outside this interval reflection coefficient /3 quickly 
decreases to zero. 

The field if can be expanded in terms of creation and annihilation operators as 



if = 



J2 ("*^-/' 



m* ikx 



lA fin -ikx 
°kJk e 



and the Green function is equal to 

G{xuh\x2,t2) = ^n{0\T^{Xl,hMx2,t2r\0)^n = j fT [t ^) Jl!'* {t>)e''^^-'---Uk = 



9{h 



to, Xi 



X2) 



The first factor here is a gauge dependent phase which must cancel in physical quantities. The 
remaining part 17(^1 — t2, xi — X2) is invariant under translations and defines correlation functions 
of gauge invariant quantities. In particular we can calculate the induced current which can be used 
to estimate the back reaction. The general formula for the current is 

{J{t)) = ( [k- A{t))\^Uk,t)\'' dk 
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As we will see, the current is dominated by the two semi-classical domains described above. Before 
the 'horizon crossing' we have 

oo 



{j{t)y'> = / dk dk^ ' = / — , ^ (6) 



A(t)<k 

where p = k — A{t) is 'physical momentum'. After horizon crossing, we have to use Keeping 
only non-oscillating terms, which are dominant, we obtain 

(j)^'^ = I dkdk^ {\a{k)\' + \mn 

J \/{k — A)'' + k, + 

k<Ait) 

Using the general relation |q;(/c)P — = 1 we get: 



J 2y/p2 + kl + m^ J 2y'p^ + kl + 

— oo — oo 

The first term in this formula combines with ([6]) and gives zero due to p ^ —p symmetry. The 
second term is really interesting. The key feature of it is that the reflection coefficient /3 depends 
on the 'comoving' momentum k and not the physical one p. As we saw, this coefficient keeps being 
constant for A{~oo) <^ k ^ ^(oo) and quickly vanishes outside this interval. In terms of p, this 
means the time-dependent cut-off A{—oo) <^ p + A{t) ^ ^(cxi). We also have a cut-off on k±, 
k± <C E. Hence, the total current is given by 



(J) = J dp^^Jdk^\l3{kj^,k)\^ = -(^A{t)-A{-oo)'j\l3\^E^-const (7) 

A(-oo)-A(t) 

In the last expression |/3| — e ^ . This result is physically transparent. It means that, as 

the time goes by, more and more k modes cross the horizon k = A{t) and begin to contribute to 

the induced current. This fact is important. It shows that the induced current is proportional 

to the vector potential and not the field strength. Together with gauge invariance this implies 

t 

highly non-local behavior. Indeed, A{t) — oo) = J dt'E{t'). Similar non-localities are well 

— oo 

known - the London equation in superconductors, the photon mass in the Schwinger Model, the 
Chern-Simons terms in the quantum Hall effect. In all these cases the gauge invariant expressions 
can't be expressed locally in terms of the the field strengths. 

The stunning feature of the above result is that the non-locality appears in the massive theory. 
This is specific for unstable vacua and can't be seen in the in/out formalism. It also implies 
the strong back reaction, since the current is growing with time. Another interpretation of this 
result is symmetry breaking. Indeed, in the constant field we obviously have the time translation 
invariance. This invariance is broken in the expression for the current due to the influence of the 
past when the field was turning on. We will return to this phenomenon later, while discussing the 
gravitational case. 

In the mean field approximation we plug the current back into Maxwell's equation similarly to 
the procedure of ^0 Considering for simplicity l+I dimensions we get 

A = J = -2\l3{A)\'^A{t) 



^In the recent paper [7] the in/in current has been calculated. The result is different from ours and we disagree 
with the method used in this paper. Another expression for the in/in current, which is consistent with our result, 
can be found in [8]. Let us also notice the paper [9] on R-N black holes which used the method very close to ours. 
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At large time the solution decays a.s E — A ^ indicating the total screening of the electric 
field. This may be related to the old result by Gribov [6 . Another analogy is the Landauer 
conductance in the mesoscopic system, which is also expressed in terms of the transmission and 
reflection coefficients. 

Another interesting point is that for massless fermions the anomaly equation reads as 

d„Ji = -E 

TT 

since axial and vector currents are related by epsilon symbol in two dimensions. This gives linear 
growth of current with time and is consistent with our result since in massless case /3 1. The 
back-reaction can be estimated from the Maxwell equation E + Ji = and is clearly significant. 
In the case of fermions (the Schwinger Model) we know that in the vacuum no electric field remains 
and electric charges are completely screened. 
We can also use the in/out Green function 

Gin/out ^ \,^out^{, \ 

a 

The sign of vacuum instability here is ImG{t\t) ^ 0. Let us notice that the matrix element 
(put\Jx\in) — simply because the in/out Green function is Lorentz invariant (module a phase 
factor). The Euclidean version of this phenomenon (with the replacement of the electric field by 
the magnetic one) is the absence of bulk currents in the quantum Hall effect. 

It is also instructive to change the gauge. If we take = Ez we get Klein-Gordon equation 

(pl + {i^-Ezf - m^)ip ^ 

As in the time-dependent gauge, we have a Schrodinger equation for inverted oscillator, but this 
time the effect of pair creation comes from the underbarrier penetration rather than from the 
overbarrier refiection. The two are related by the analytic continuation. In this gauge the energy 
UJ = Ez -f \/ + TTi^ is conserved but non-positive which allows particle production. This 
gauge has de Sitter counterpart. The tunneling above is analogous to the tunneling in the Painleve 
coordinates |10| . 

4 Expanding Universe (free fields) 

When we look at the de Sitter space, we find that there are striking similarities with the electric 
case. Let us consider what happens when the curvature of dS space is adiabatically switched on. 
In this setting we have two quite different problems - expanding and contracting universes. The 
arrow of time is set up by defining the infinite past as a Minkowski space in which our field is in 
the ground state and solutions to the wave equation are chosen to be the Jost functions. Let us 
begin with the expanding Universe. Analogously to the electric case we will assume that the FRW 
metric 

ds^ = a{tfdx^ - dt^ 

is such that ^ = H{^)^ time T is supposed to be large, and Hifi) = 1, while H{±oo) = 0. A 
representative example of such a metric is 

a{t) = e'^t^^h^ 

H{t) = cos}i(-^Y • ^® convenient to rescale the standard scalar field Lp by defining ip = a~^(j). 
The Klein-Gordon equation takes the form 
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with r{t) = ^^^(f )2 + f f . As before, the 'in' solution is defined by 

^ —iuj.t 



e"^". * (8) 



as t — > — oo with cj^. = Im^ + ^^_.^yi j ■ Its quasiclassical expression is given by the formula 



where iOk{t) = w'm^ — r(t) + This WKB expression is applicable if 



d 



(1- 

dt \ U!k 



< 1 



If we assume that H ^ - ^ m and H is small, we see that WKB breaks down when the given 
mode crosses the horizon, k ^ ma{t). Before that we had k ^ ma{t) and A ^ 1. Long after that 
we reach the semi-classical regime again, but with two exponentials as in ([3]). Let us consider the 
time evolution of the quantity {in\(p{t)'^\in). We have: 

{in\^{tf\in} = I d'^k |<y9„(t,fc)|2 



splitting the integral as before into the regions |fc| S> ma{t) and |fc| <C ma{t) we get 



2uJk{t) ' ' J 2uk{t) 

|fe|>ma(t) \k\-^ma(t) 
\k\<^7na{t) 



(9) 



The reflection amplitude /3(fc) is fc-independent in a certain interval, just as it was in the electric 
case. The reason is that de Sitter wave equation is invariant under fc — )■ Afc and t — >■ t + log A 
(which is one of the dS isometrics). However this amplitude quickly vanishes when k is such 
that the horizon crossing happens outside the de Sitter stage. Namely, if tk is determined from 
the equation k = ma{tk), the constant reflection occurs for \tk\ T. If we introduce the cut- 
offs defined by ^^""^ = kmin^^ = ^ a-nd ^^"^j = kmax^~'^ — we have reflection only if 
kmin <C fc ^ kmax- We sec that the contribution of the second term in ([9]), which represents the 
created particles, is small in the expanding case. Due to the infrared convergence of the integral 
we obtain 

{ip{t)')<^'^ ^ |/3|2m'^-i (10) 

This formula has a clear physical interpretation. By the moment t we excite the modes with 
< mait) and the average excitation number is n ~ The created particles are non- 

relativistic due to the upper boundary on k. Let us stress that there is no dilution of the created 
particles in the sense that their physical (not comoving) density remains constant in time, however 
their main contribution is just a renormalization of the cosmological constant which is unobservable. 

The key difference from the electric case is the absence of the dynamical symmetry breaking, 
which we define as a long-term memory. By this we mean the following. As we already noticed, 
the current in the electric case depends on the time passed from the first appearance of the field. 
This effect is a dynamical counterpart of the usual spontaneous symmetry breaking. In the latter 
case, the magnetic field at the boundary induces magnetic moment in the bulk, if we talk about 
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Ising model for example. In our case the role of the boundary is played by the infinite past. The 
expression (jlOp does not depend on time. Hence, there is no dynamical breaking of de Sitter 
symmetry in this case. Life becomes more interesting if we switch on interaction or consider 
contracting universe. 

We could calculate things in the regime of the constant curvature and get the right results. In 
this case 

with T — e~* and 

{^nW{tf\^n) ^ j k \h[]^ {kT)\^ = jd'p\H\^{p)? = const 

The UV divergence in this integral is the same as in the flat space and the time independence 
in this formula is just the result of the de Sitter symmetry. The back reaction is thus small and 
uninteresting. Really non-trivial things begin to happen when we either include interactions or 
consider contracting universe. We start with the latter. 



5 Contracting Universe (free fields) 

Let us repeat the above calculations in the case of contracting Universal. At the first glance it 
may seem that, since de Sitter space is time-symmetric, expansion and contraction can't lead to 
different results. However, as was stated above there is an arrow of time in our problem. We 
defined the past by the condition that our field is in the Minkowsky vacuum state. Generally 
speaking, in the future we should expect complicated excited state. In this setting contraction is 
very different from expansion. We can once again take 



a{t) 



-Tta 



The modes with k > ma (—00) 



will always stay in the WKB regime, since a{t) will be 



decreasing. On the other hand, the modes with ma{oo) <C A: ^ ■ma{— 00) will cross the horizon 
at some time, k « ma(tk). If we once again define the 'in' modes, ipin{k,t) by the condition ([8|, 
we find that for k ^ ma{t) the horizon crossing (WKB breaking) has not occurred yet (remember 
that a{t) is decreasing) and hence we have a single exponential ([8]). 

For ma(t) <^ k <^ ma{— 00) the horizon crossing is already in the past and we have two 
exponentials with the coefficients a and /3, |q;(A:)P — |/3(fc)p = 1. For fc ^ TOa(—cx)), the horizon 
crossing has never occurred and /3 — ^ 0. As in the previous section we get 



{in\Lp{t)'^\in) — a{t) 



d'^k 
2oJk{t) 

\k\<^7na{t) . \k\^ma( — 00) 

r d'^k 



2L0k{t) J 2ujk{t) 

rna{t)<^\k\<^rna{-~ 00) 



, , f d'^k 
^ ' J 2uJkit) 

ma{t) <^\k\<^ma( — 00) 



(11) 



Collecting different terms we get 



d'^k 
2ujk 



|fc|<Aa(t) 

const ■ A"*"^ + 



d'^k 
2uJkit) 

ma{t)<\k\<ma{ — 00) 

a{t) J 



(12) 



We considerd the case of contracting universe following the advice of V. Mukhanov. 
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The first term in this formula is just the same UV divergent term as in the Minkowsky space. 
The heart of the matter is the second term which displays the symmetry breaking through the 
long-term memory (dependence on a(— oo)). However, the memory can't be too long, since we 
have a standard UV cut-off at the Planck mass. Because of it, the above formulae are valid if 
p = k/a{t) < Mpi and therefore a{—oo)/a{t) < Mpi/m. 

Let us sum up the above discussion. In the expanding universe the contribution from the 
created particles comes from the region ma(— oo) <C fc <C ma{t). No long term memory is present 

and time-dependent back reaction is small, of the order of (j^^'^j^^ ■ Created particles are 
non-relativistic due to the red shift. 

In the case of contracting universe particles come from ma{t) < < min (rna{—oo), Mpia{t)'j . 

They are ultra- relativistic and their contribution is of the order ^ "^^^^^-' ^ — > oo. All these 
conclusions are correct only for non-interacting particles. 

It is also possible to calculate the energy-momentum tensor. We have 



In the contracting case the order of magnitude of this quantity is defined by the integral: 



This corresponds to the ultra- relativistic particles with the equation of state p — ^e. In the 
expanding case the contribution to Tqq comes from a small number of created non-relativistic 
particles. In both cases there are no reasons to believe that created particles are in thermal 
equilibrium. Let us also stress that the above formula represents a non-local contribution to Too 
similar to ([7]) . In contrast with this formula, the local contributions should depend on the quantities 
taken at the time t only. 

6 Secular interactions and the leading logarithms, Poincare 
patch 

In this section we discuss a very peculiar property of the de Sitter space. Namely, it turns out 
that the interactions of the massive particles generate infrared corrections. We start with the 
second order of perturbation theory in the case of Xip'^ interactions (which we choose to simplify 
notations; the phenomenon we are after is general and has nothing to do with the naive lack of the 
ground state of the above interaction). We first calculate the correction to the Green's function 
G{q,T) = {in\ip(q,T)ip{—q,T)\in) where q is a, comoving momentum in the Poincare patch and r 
is a conformal time. Our goal is to show that if the physical momentum p = qr <^ fj,, there are 
corrections of the order (A^ log ^)" where /z is the particle mass; notice also that these logarithms 
are the powers of the physical time t — — logr. 

We are interested in the loop corrections to the one-point function {ip{t)'^). The magnitude of 
this quantity determines the strength of the backreaction. To find it we have to use the Schwinger- 
Keldysh perturbation theory. These methods are well known and we will add a few explanations 
to fix notations. Let us suppress first the momentum dependence and expand if = f*a + fa~^ , 
where f(t) are the "in" modes and a is an annihilation operator. The relevant one-loop diagram 





ma<k<'ma{ — oo) 
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is shown at Fig. 3. It's contribution to G{q, r) ~ {in\ ip{q, r) ip{—q, t) \in) is given by 

t ^ 

G{q,T) = -x^ r.itf J dhdh f,{h)f,{h) J ^ h{t<)m>)h+s<)rk+,{t>) - cc. + 

— oo 

t 

+ 2.A2|/,(t)p j dhdt, f,it,)f;{t,) j f^{t,)ft{h)h+,{ii)n+,{h) 

— oo 

(13) 

In the first line we have the contribution of the (+/+) and (— /— ) diagrams (the signs refer to the 
points ti^2 of the physical time, or ri.2 of conformal time at the diagram in Fig. 3), while in the 
second line we have (+/— ) and (— /+) diagrams. 




Figure 3: One- loop diagram responsible for infrared logarithms in Poincare patch. 

We choose the "in" wave function to be a Hankel function 

fk{t) ^ rihikr) = constT'^/^H[^^{kT) 

where the normalization is fixed by the condition h{x) (2a;)~^e" as a; —J' oo. With this 
normalization, the asymptotic behavior at a; —> is given by 

h{x) A{^l)x'^' + A{~^l)x~'^' (14) 

where As are some concrete functions which we discuss later. 

As we will show in a moment, there are infrared logarithmic corrections to G{q,T) = r'^giqr) 
when qr ^ /i. In this regime we can use asymptotic expressions to get, A = log (^) 

gix) = A{ii)A*{~fi)rCX,t^)x^'^ + Ai-t,)A*{fi)T*Cx,fi)x-^'^ + {\Aif,)\' + \A{-ti)\') C(A,/i) 

when interaction is off (A = 0), coefficients F*^"-' — C^^^ ~ 1. Our goal is to find these quantities 
at non-zero A. We start with the interference term C. 

In order to obtain the contribution to giffr) we have to integrate the diagrams of Fig. 3 over 
the momentum k and the time variables t\ and ^2 • The logarithmic contribution comes from the 
domain Ti 2 ^ [ijk and ^/t ^ k ^ q. In this domain we get the contribution from the first term 
in in the form 

//■oo /-oo 
d'^k j dTlJ dT2{TlT2f/^~^h{qTl)h{qT2)h*{kTlfh{kT2f-C.C. 
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Taking the limit g — > and interchanging 1 and 2 in the complex conjugate term gives 

<?M/ = ^^^c,(M) = c,iog(^) 

Here the coefficient is given by 

Ci = -4A2|A(M)A(-/.)|2(|g(^)|2 + |g(_;,)|2) 

POO 

g(p) = dx e{x) 

Jq 

The second term is treated analogously. It has the form 

g{qT)ii = 2\^h*{qT)h{qT) j d'^k dTidT2{TiT2Y'^-^h* {qTl)h{qT2)h* {kTlfh{kT2f (15) 

Integration gives another logarithm. Summing these contributions finally gives for the interference 
term 

{vD = 9{qT)i+g{qT)u = 2.(B(^)-i?(-M))-(s(M)l5(M)l'-S(-/i)l5(-M)l')-A2log(^) 
where 

Bill) = \A{y)\^ = e"^ . , 
^ ' ' 4/i smh(7r/i) 

The first multiple here is a Wronskian of the eigenmodes. The second one turns out to be equal 
to zero. To see this, note that the functions h{x) satisfy h{x)* — e^^ h{e^'^x) which implies the 
following relation for g{fi): 

\g{pr = e-2-'%(-M)P 

The physical meaning of this equality is detailed balance relation with Gibbons-Hawking temper- 
ature for de Sitter space. Combining this with the similar property for A{fi), we conclude that the 
one-loop contribution to the coefficient in front of the logarithmic divergence in the interference 
term is equal to zero C^^''(A, /i) = 0. 

The next step is to calculate F. Imaginary part of this quantity determines the renormalization 
of mass ^, which we are not interested in at the moment. The real part is responsible for the 
imaginary contribution to /i, which is related to the decay rate of the particle. Using similar 
trickflfl to those used above we find 

i?e(rW) = A2(B(M)-i?(-M))(|5(M)P-|5(-M)P)log(^) 

This quantity is non-zero and negative. 

The above calculation refers to the IR properties of the two-point function. In the case of 
Poincare patch there is no IR contribution to the one-point quantities, as can be seen from the 
conformal diagram at Fig. 4. The Poincare patch is shown here by the gray area. Interactions 
contributing to the one-point function must be located inside the past light cone due to causality. 
Therefore we have to consider only the intersection of the light-cone with the gray area defining 
Poincare patch. Thus infrared effects in Poincare patch can not have dramatic consequences 

''It is convenient to rescale k from the integrals over ri,2 and note that 

oo oo 

Y = jdxj dy {xy)i-^[{^Y^ + {"^y^^Yiyfh'ixf = \{\9{tA? + \g{-^A?) + iA 

X 

where A is some real number, contributing to renormalization of ^ only. 
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Figure 4: Conformal diagram. Poincare patch is shown by the gray area. Solid black line represents the past light 
cone of the observer. The intersection of this cone with Poincare patch touches past infinity only at one point. 



because the past iufinity is represented only by one point. In the complete de Sitter space the 
situation is quite different and is discussed in the next section. 

Although infrared corrections do not appear in the 1-point function {(p{t)'^), they contribute 
to the two point function {ip(l)ip(2)) . To illustrate this consider the limit when ri = T2 = t and 
x'^ = {xi — 2:2)^ ^ T^. This corresponds to z ^ —00 (|23p . The bare Green's function in this limit 
is given by 



Goiz,fj.) = 



1 



The exact Green's function is equal tc0 



G{z) 



1 + ^(B{^,)-B{-^,)){\g{^i)\^ 



1 



A2 
4/i 



(1 - e-^-A') |.g(-M) P log(-z) Go(z, + 6f,) 



We see that besides the infrared renormalization of mass, which we ignore in the present paper, 
the bare Green's function is multiplied by the function of log(— z). Thus, even in Poincare patch, 
infrared corrections do appear when the two points are separated by a large geodesic distance. 
It would be interesting to understand the consequences of this result for the inflationary models 
in Poincare patch. 



7 Secular interactions and leading logarithms, complete dS 
space 

In order to describe the global dS space, we use the standard metric ds^ ~ dt^ — cosh^ t{d£ldY ■ The 
eigenmodes for the Bunch-Davies vacuum are inherited from the sphere. To simplify notations we 
write them for d = 1: 

fg{t)^ PZl^^Ji sinht) 

where q is an integer. These modes are selected by the condition that they are regular when 
continued to the southern hemisphere {t = —M; d > 0). 

■'To derive this formula we can make a Fourier transform 
i£ 

1^ dq-T[ A(^,) A' i-f,)r {qrf^^ + A* (/.) T* (gr)-^^^ + {\A{^)\^ + \A{~t,)\^) C ] e'"- 

and retain only terms of the order log(— z) while neglecting the terms of the order . 
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The logarithmic divergences appear when j^j 3> 1 and |i| — >■ oo. In this regions the Legendre 
functions can be replaced by the Bessel functions. We have: 

•^'^ %^oo |r'^/2/j*(gf), T = e+\ t^-oo. 

As it should be, this is exactly the doubled Poincare patch. 

Let us use these modes to calculate perturbative corrections to {ip'^{n)), assuming that the 
interaction begins adiabatically in the far past, with t = £ — > 0, while the "observer" sits in the 
future at fixed r. The most important contribution comes from the H — term in the Fig. 5. We 




Figure 5: Relevant diagram, leading to IR divergence, in complete dS space. 

have 

oo 

{^'in))('^ = ^^\h{qr)\'J ^ {h* (qn) Hgr,)) ■ (nr,)^ -a.inj,) (16) 

e 

where 

/d'^k _ _ _ _ 

j^-^ h*{kTi) h{kT2)-h*{\k - q\Ti) h{\k - q\T2). 

We consider here only the dominant contribution, when ti, t2 are both in the far past. If fc 3> 
we get the following scaling property: 



'q 

The integral becomes: 



oo 



€ 

= const ■ \^ -T-^j d-^q\h{qT)\^ \og[^^ 
The UV divergence at large q must be cut-off by the condition qT < Mpi . Thus we get the result 
= const ■\H4-{Hog[-±-^''-) = const ■ \^ {^\n))^^) ■ \og[-^^ . (17) 
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This formula is valid if: 

A* 

IVIpi 

which means that the time T during which the interaction was on, satisfies 

1 ,T, 1 /"Mpi 

T= - log - » 77 log — 

(where we reinstated the Hubble constant). 

In the Schwinger - Keldysh language we accounted for the (+/— ) self-energy part. There 
are, of course other insertions, (+/+) and (—/—), also generating secular logarithms. However, 
they are proportional to J d'^q h?{qT) and its conjugate. This integral is UV convergent due to 
the oscillations of h{qT). Hence there are no UV/IR mixing in these terms and their secular 
contribution, while non-zero, does not contain Mpi, unlike (I17[) . 

It is also instructive to write the above correction in the covariant form. Various Schwinger- 
Keldysh propagators are expressed in terms of the different boundary values of a single analytic 

Wightman function, g{n ■ n'), e.g. (?++ = g{n ■ n' — «0), = g{n ■ n' + iesgn{nQ — tiq)) etc. 

The function g{z) is real for 2; < 1 which corresponds to the space-like separations. By combining 
terms in the Schwinger- Keldysh diagrams it is easy to get {ip'^) = {lp^)i + {^p^)ii wher^l 



((/S^)/ ~ j dnidn2 lTn(^g{nni)j Im(^g{nn2)j ({ J(ni), J(n2)}) 
{ip'^)n -2i I dnidn2 0{nio - n2o) Im{g{nni)^ Re{g{nn2)^ {[J {ni) , J {n2)\) 



where J ^P' (n) . The first term represents the contribution of the real particles created from the 
vacuum, while the second term comes from the virtual particles. The logarithmic divergence arises 
from the domain where nio, n20 ^ —00, while (nin2) ~ 1. The mathematical origin of the UV/IR 
mixing lies in the fact that the first integrand contains terms g(z + iO)g{z' — iO) which become 
singular on the light cone while infrared divergent in the infinite past. This phenomenon never 
happens in the Minkowski space. 

In higher orders there are higher powers of the logarithms. Their summation requires the 
kinetic equation and will be discussed elsewhere. 



8 Conclusions 

The physical interpretation of the above estimates is the following. We are considering a complete 
dS space. All points of this space are equivalent, so that the statements that at a given point we 
have expansion or contraction are meaningless. However, if we fix the position of the observer, one 
can define domains, such that the signal sent from them will be either red shifted or blue shifted. 
The essence of the formula (|17l) can be grasped from the Fig. 5. We integrate the interaction over 
the faraway past region. The size of the loop determines the interaction scale ^ 1/m, which is a 
large quantity. While the signal from the interaction region propagates along the geodesies to the 
observer, sitting at the point r, it is blue-shifted to the Planck scale ^ 1/Mpi. As a result we get 

^The easiest way to derive these formulas is to use the definition of (anti)chronological products 
T{.hJ2) = i{Jl,J2} + 'rhjuJi], f(JiJ2) = i{Ji, J2} - -rilJi, J2], J1J2 = i{Ji,J2} + hh,J2], 

where u = sign{nio — n2o) and then use the symmetry of the measure w.r.t. intrchange of 1 2 to reduce the 
integration domain to the region nio > n2o- AH the terms containing anticommutator of currents are collected into 
{ip^)l, all the terms containing commutator are in {ip^)ij. 
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a very curious UV/IR mixing. In the fiat space we expect tliat UV and IR divergencies contribute 
to the physical quantities independently - we do not expect the terms, like (I17[) . which are both 
UV and IR divergent at the same time. This is a specific feature of the curved space. 

The £ dependence of the physical matrix element discussed above indicates a breakdown of the 
dS symmetry; as always, spontaneous symmetry breaking manifests itself through the sensitive 
dependence on the boundary conditions. The logarithms will be present even for a patch of the dS 
space, provided that it is "large", that is the past cone of the observer intersects a decent portion 
of the past infinity. As we saw from Fig. 4, this is not the case for the Poincare patch; for it "the 
world is not enough". 

Let us explain our motives for using the global dS space, while in the inflationary theories 
only a small part of it is usually present. Our goal is to resolve the puzzle of the cosmological 
constant by infrared means. We start with the Einstein action with the cosmological constant 
present. The standard procedure in field theory is to assume first that we can neglect quantum 
corrections at large distances, find a classical solution and then evaluate the corrections. It is this 
procedure which allows us to use classical Einstein or Navier - Stokes and forbids the similar use 
of the Yang - Mills equations (due to asymptotic freedom) and sometimes the diffusion equation 
(due to Anderson's localization). 

In such a setting we must consider the global dS space as a first step. If a starting point were 
incomplete space, we would end up with the unitarity problem, since the particles can disappear 
from the space. Of course it is possible to have a Poincare patch glued to the Minkowski one in 
such a way that the result is geodesically complete. However this space will not be a solution of the 
Einstein equations with the cosmological constant only. It is also possible to modify the Einstein 
action so that we have a different background without IR divergences. This looks ambiguous and 
is far from our goal, which is to tame infrared divergencies. We should remember that to solve the 
A-problem one must be searching for the infrared effects and not running from them. IR divergence 
is not a problem but an opportunity. 

Another question is related to the choice of the Bunch-Davies vacuum in the above calculation. 
What is the reason for this (apart from the tradition)? It seems that the right starting point 
should be the state with the longest life time. We haven't proved that this is the case, but various 
estimates make us believe that the Bunch-Davies vacuum is the most stable one. In the appendix 
we present the propagators for the different possible ground states. It should not be difficult to 
extend our analysis to other vacua. 

Finally, there is a number of valuable papers |15| intersecting with our work, but it seems that 
our approach brought some new and unusual results. 
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Appendix A. Oscillator, relative probabilities 

Here we briefly discuss some ideas mentioned in the main text using the simplest model - quantum 
mechanical particle. Take an oscillator with variable frequency 

{d^ + + U{t))(p = 
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Let us introduce Jost functions 



fout{t) -> 



/2m 



/2m 



t — ^ —oo 
t -)■ +00 



(18) 



As well known in scattering theory 

fUt) = afoudt) + PfLtit) 

|ap — |/3p = 1. Let us find the vacuum decaying amplitude. We define the (in) and (out) vacua 
in a usual way: 

and a\0)in — b\0)out — 0, [a, a^] = [6, = 1. The in/out Green's function 

1 1 



G 

— out 



fin{t<)fout{'t>) 



out{0\0)in a 

This Green's function is the one satisfying the composition principle [T] 



dw? 



i I G{h,t)G{tM)dt 



This equation allows to represent G in terms of the Feynman's sum over paths. The amplitude to 
produce 2n-particles is given by: 



C2n 



Ao-^2n = out(0|^=|0)j„ = o«i(2n|0)i 

v2n! 



Let us express it in terms of the Green functions: 



G(tl, t2n) — 



out (0|r^(ii)...(^(t2„)|0) in out \^ I " in \ 



^tj —¥oo 



t{m^n (2m)'^ 



out 

(0|0) 

in 

The Wick theorem on another hand gives 

G{ti,...t2n) = G{ti,t2)-G{t2n~i,t2n) +perm. 
Total number of permutations is (2?! — l)(2n — 3)... = (2n — 1)!! = From here we derive 



(2n)!//3\" 1 



2"n 



(we pick up the terms containing e from the products of the Green's functions). Hence the 

probability to produce 2n particles is 



W' 



\oiLt 



(2n)! 



Wn 



The normalization condition gives 

00 00 

I = Wo + J2W2n ^ Wo{l + J2 



o„t(0|0),„|' 
~" W2n^ 



(n!)24" 



n=l 



W, 







(2n)! 



(n!)24» 



2n 
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Thus Wq = yl — I^P = It is interesting to note that the Green functions give relative 

probabihties, Po„ = When the vacuum is unstable Wq < 1 and hence J2 Pon > 1- This may 

n 

open the way to interpretation of the non-unitary field theories - they describe unstable vacua. 
For example in non-unitary CFT there are well classified operators with negative norms. The 
'probabilities' extracted from the Green functions satisfy ^±P„ = 1, > 1 and P„ can be 

n n 

interpreted as relative probabilities. 

Another interesting relation is Wo = e^^, ^ J Im G{t,t) dt = Im(^^ ^ f^utiifdtj . We 
see that Im G{t, t) signals vacuum instability. The typical back-reaction of produced particles 
on the field U{t) is characterized by the 'current' J = in{0\v{t)'^\^)in- We see that 

J{t) = \fin{t)\' 

and (omitting rapidly oscillating terms) 

J(oo)-J(-oo) = i^(|a|2 + |^|2-l) = 

zm m. 

Thus, the back-reaction docs not become large with time. We can also evaluate the average number 
of the produced excitations n = X^W2n(2n) = However it is to be remembered that this 

number fluctuates: (n — n)^ n. 



Appendix B. The Green's functions of the de Sitter space. 

The propagator is a solution to the inhomogeneous Klein-Gordon equation with a delta-function 
source. Such solution is not unique and defined modulo a solution to homogeneous equation. 
In the flat Minkowski space there is a well established prescription of doing integration in the 
complex momentum plane which fixes this ambiguity and uniquely defines the propagators. This 
prescription however is not straightforwardly generalized to curved space which we are working 
with, hence we need a different setup. A convenient method is to write the propagator in the form 

Where (^1,2 (i) are two linearly independent solutions to the homogeneous equation, W[<pi,(/?2] 
is their Wronskian and f<> is the smallest and largest of the times ti_2. The above mentioned 
ambiguity is reflected here in the different possible choices for the solutions to the homogeneous 
equation <^i,2(*)- This choice depends on the problem we want to solve. In case of in-out propagator 




Figure 6: On the left - geometry of expanding Universe, on the right - geometry of complete de Sitter space, 
we should choose v'i,2(i) to be Jost functions of the scattering problem, i.e. ipi to be a plane wave 
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at t — ^ — oo and (/32 to be a plane wave at t ^> +00. In case of in-in propagator both ipi and ip2 
should be plane waves at t — > —00. 

Another complication comes from the fact that a de Sitter patch can be embedded into flat 
space in numerous ways. Below we illustrate this by considering two possibilities: initially flat 
space adiabatically starts to expand; and complete de Sitter space glued with two cylinders on 
both plus and minus time infinities (see Fig. 6). The Green's functions for these geometries turn 
out to be very different. 



9.1 First geometry, Poincare patch 

We start with the first geometry, which can be conveniently described by the Poincare patch 
ds^ — (dr^ — dx'^) and for simplicity we work in d = 1 dimensions. The wave equation for 
a mode with comoving momentum p can be written as 



with u = j/i = (ma)^ — i. It's solution can be written in terms of Bessel functions. In the 
in-region the Jost functions are Hankel functions, while in the out region they are Bessel functions. 
Thus, according to (|19p we obtain 

00 

G^p"/°"*(1,2) = J dp cos[pix2-xi)]V^ Hl'\pT^)J,{pr^) (20) 


for in-out (Feynman) propagator and 

00 

G*p"/"(1, 2) - J dp cos[p(x2 - X,)] ^ Hl'\pT^)Hi^\pT^) (21) 


for in-in propagator. Integrals w.r.t. momentum can be carried out explicitly by using 6.672.4, 
6.672.3 of [Ti]. The resuh is 

Gp/°"*(1,2) = Q,_iiz + iO) (22) 

G™/'"(1,2) = -^V^_^i-z-^0) 
cos(7ri^) 2 

The first expression was suggested in [IJ, the second one is the Bunch-Davies propagator. Both 
are expressed in terms of geodesic distance 

T^+rl - ix2-xiy 
2riT2 

To make the integrals (P0)l . pT|) convergent we need to shift r> into complex plane. This shift 
together with > determines iO prescription in the arguments of Legendre functions. 
Another nice representation of in-out propagator can be obtained from 6.669.3-4 of |14) 

00 00 

G™/°"*(1,2) = J dp cos[px2i]V^e--^ J ^/''^'' + -^^"{^^v) ^^^^ 
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This expression are analogous to Schwinger proper time representation of the Feynman propagator 
in electric field. It's path integral derivation can be found in |^13j. It is also convenient for the 
calculation of the imaginary part of in-out propagator at coincident points, which is related to the 
imaginary part of effective action. Taking the limit of coincident points in (j24p and integrating 
over k by using 6.611.1 we get up to inessential constant 



Im 



Gp/°"* (1,1)' 



+ 00 ^ 

^ Im— f ^^e^'^^ = -y(-l)"e-2™M ^ 
2tt / sinhs ^ 



-27TfJ, 



1 + 6-2^^^ 



Non- vanishing of this quantity signals instability of the vacuum w.r.t. creation of particles, similarly 
to Schwinger mechanism in constant electric field. 



9.2 Geometry of the Complete de Sitter space 

The second geometry can be described by a metric 

ds'^ = dt'^ - co&ht^dip^ 

with compact coordinate S [0, 2-k]. The wave equation for a mode with integer momentum p can 
be written as 



+ tanht0 + m (/) + 



P 



cosh t'^ 







The general solution of this equation i^ 

1 



CiP±^i(±tanhi) + C2Q*^i(±tanht) 

'P 2 P 2 



Vcoshi 

Choosing the Jost functions and plugging them into (fT9)) we obtain 

00 

-iiii/in 



G 



dS 



2-\/coshii cosht2 



^ Ep cos {pip) tanh t > ) Pp_ j. (- tanh t< ) 



p=0 



a 



■in I out 
dS 



1 



Vcoshii cosht2 o 



00 ^ ^ 

Y^epr{v-p+-)V{v+p+-) cos(p^)p-^ (tanht>)p-^, (- tanhi<) 



Ep cos{p(p) Pp^ 1 (tanh ty ) Q'^_ 1 (tanh t< ) 



cos(7r:^)-\/coshti coshi2 

where gp = 1 for p = and £p = 2 otherwise. It is possible to sum up the p-series and express the 
result in terms of Lorentz invariant quantities (geodesic distance z and a =^ Sign[no{l) + no(2)]). 



G 



in I in 
dS 



1 r 



Q,i{-z-ie) - Q,i{-z + ie) 



4j cos(7ri/) 



(a + 1) V,_dz + ie) - V,_.{z-ie) 



We would like to emphasize that this expression is Lorentz invariant, since functions Vj^_i{z ±ie) 
have a cut only for z < — 1, but cr is Lorentz invariant quantity for z < — 1. Also, this function 



^In the following we adopt the notation P'i{x) and Q!i(x) for associated Legendre functions on the cut when the 
argument is —1 < < 1. These functions are defined by |14j 8.702-8.705. For Legendre functions in the complex 
plane we use symbols Vj^lz) and Q{!^{z). 
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vanishes within the hght-cone of the past of the antipocQ- This area corresponds to z < — 1, 
a = —1, see Fig. 7. 




Note, that this answer is different from in-out propagator in Poincare patch (P^ . c.f. also ^Jj. It 
also has a non- vanishing imaginary part at coincident points. 
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